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O . Abstract 

^■"^ ■ By means of a novel variational approach we study ergodic properties of a model of a 

multi lane traffic flow, considered as a (deterministic) wandering of interacting particles on 
an infinite lattice. For a class of initial configurations of particles (roughly speaking satis- 
fying the Law of Large Numbers) the complete description of their limit (in time) behavior 
is obtained, as well as estimates of the transient period, in this period the main object of 
interest is the dynamics of 'traffic jams', which is rigorously defined and studied. It is shown 
that the dynamical system under consideration is chaotic in a sense that its topological en- 
tropy (calculated explicitly) is positive. Statistical quantities describing limit configurations 
C*^ ■ are obtained as well 

o 
o 
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> ! 1 Introduction 
v< 

Despite a self evident practical importance of the analysis of traffic flows and a relatively long 
history of attempts of their scientific treatment (going back to the fifties) only recently (in the 
end of nineties) reasonable mathematical models of traffic flows and method to study them 
were obtained. Previous attempts were based on ideas borrowed from such classical fields of 
physics as mechanics and hydrodynamics. Not going into details of a qualitative and quantitative 
comparison of the hydrodynamic type models with practice (which one can find, for example, 
in a recent review and references therein), we consider the following practical observation. 
It turns out that going by foot in a slowly moving crowd it is faster to go against the "flow" 
than in the same direction as other people go. A mathematical model describing this effect 
in the case of the one lane traffic was introduced in [gj. A standard probabilistic model of a 
diffusion of particles against/along the flow clearly contradicts to this observation, which very 
likely indicates a very special (nonrandom) intrinsic structure of the flow in this case. The main 
aim of the present paper is to study how this structure emerges from arbitrary (random) initial 
configurations in a simple model of the multi lane traffic flow. 

A recent progress in the analysis of traffic flows was due to the introduction of discrete (in 
time and in space) cellular automata models of the one lane traffic flow in []i~4| , |i~3|| and later 
studied by many authors (see Q| for review and further references). Various approaches starting 
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from the mean-field approximation |TT|] to combinatorial techniques and statistical mechanics 
methods [llj were used in their analysis. All these models were based on the idea to describe the 
dynamics in terms of deterministic or random cellular automata (see results about stochastic 
models in |], [7|, ||, ||]) and to a large extent were studied by means of numerical simulation 
(especially because of low computational cost of the numerical realization of cellular automata 
rules, which made it possible to realize large-scale real-time simulations of urban traffic p0|]). 

Roughly speaking the one lane road in these models is associated to a finite one-dimensional 
integer lattice of size N with periodic boundary conditions and each position on the lattice is 
either occupied by a particle (represented a vehicle), or empty. On the next time step each 
particle remains on its place if the next position is occupied, and moves forward by one place 
otherwise. In |]l4| , |i~0| , Q it was shown (mainly numerically and by some physical argument) that 
limit (as time going to infinity) behavior of the dynamical system under configuration depends 
only on the density of particles in the initial configuration. This result was generalized for the 
case of the dynamics on the infinite lattice and proved mathematically in [g], where a novel 
variational approach was introduced. 

Despite various generalizations the one lane restriction of these models was crucial, for ex- 
ample, in an attempt to study a multi lane model in Q, satisfying standard traffic rules, no 
mathematically interesting phenomena were found. Only two years ago in [15, 16] a first non- 
trivial multi lane generalization was introduced for the case of a motion on a finite lattice with 
periodic boundary conditions, based on a so called ultradiscrete limit of the well known Burgers 
equation. 

In the present paper is we study ergodic properties of this model. As we shall show this 
analysis boils down to the study of the dynamics of 'traffic jams' (see the rigorous definition in 
Section |I]), which mainly depends of the density of particles in the initial configuration. 

One of the main quantity of interest in traffic models - the average velocity of cars and its 
dependence on the density of cars p (called the fundamental diagram) is typically studied in 
the steady state. From our results it follows that in the multi lane model that we consider the 
average velocity in the steady state is equal to max{l, K/p — 1}, which immediately reminds 
the similar result known in the one lane case. 

The paper is organized as follows. In Section ^ we describe the model in detail and introduce 
the basic notation including important notions of dual configurations and maps. In Section ||] 
we introduce the space of regular (statistically defined) configurations, show that this space is 
invariant with respect to dynamics and formulate the main result of the paper - Theorem 3.1 . 
Qualitatively this result means that in a steady state any configuration either consists of free 
(moving independently) particles, or this property holds for all empty places on the lattice. In 
terms of mentioned above variational principle this can be formulated that the total number 
of free particles between two fixed ones can only grow in time. The proof of this result in the 
next Section || is based on the detailed analysis of the dynamics of traffic jams. It is worth 
notice that in distinction to the one lane case the formal description of the traffic jam is rather 
nontrivial and some individual particles in it can move still representing obstacles to the motion 
of other particles. Section |] is dedicated to the proof of chaoticity of this model: we explicitly 
calculate its topological entropy and show that it is strictly positive. In the last section || we 
derive statistical quantities describing typical limit configurations. 

We tried to define rigorously all important objects that we consider in the text, however 
of course we were unable to introduce all standard mathematical definitions. The reader can 
find exact definitions and further references related to dynamical systems (especially acting on 
discrete phase spaces), for example, in [l], |i^] . 



2 Multi lane traffic flow models: dynamics in space of configu- 
rations 

The model corresponds to the highway traffic flow on a road with K lanes. Let := 

{0, 1, . . . , K} be an infinite lattice, positions on which we call (lattice) sites. For a sequence 
X G Xq- and x G 5Z 1 by X(x) we denote the x-th element of this sequence. Consider a map 
T : Xq — > Xq", defined by the relation 

TX(x) = X(x) + min{X(x - 1), if - - min{X(x), K - X(x + 1)}. (2.1) 

Remark. In the above map were introduced for the case of the finite lattice with 

periodic boundary conditions. Observe that a finite lattice of arbitrary size N < oo with 
periodic boundary conditions is a particular case of 7L 1 lattice considered in our paper restricted 
to only iV-periodic configurations. The paper [|l(| claims an estimate of the transient period in 
the 2A r -periodic case as N. Since the construction in this paper sensitively depends on the length 
of the period it cannot be extended to the case of unbounded lattices with general nonperiodic 
initial configurations that we consider. It is of interest that 2N + 1 periodic initial configurations 
lead (as we shall show) to much worse estimate of the transient period 2N. 

The above map can be described in a different way in terms of configurations of particles. 
Let us introduce several definitions. 

A collection of particles 3 on the lattice 2Z 1 we shall call ordered if there is a function (called 
index function) I : H — > 7L l such that for any two particles £, £' G S, £ ^ £' the corresponding 
indices satisfy the inequality /(£) ^ -/"(£') and if these particles are located at sites |£| < on 
the lattice then /(£) < /(£'), where |£| stands for the location of the particle £ on the lattice. 

To a configuration X G Xq' we associate an ordered collection (finite or infinite) of particles 
on the lattice 7L l containing not more than K particles at each site, such that X(x) for x G 7L l 
means the number of particles located at the site x. Then the set of all possible ordered 
configurations of particles containing not more than K particles at each site forms the phase 
space Xq := {0, 1, ... , K} of the system under configuration. 

For a given positive integer K and a given configuration X G Xq" the action of the map T 
can be described as follows. For each site x G 7L 1 (independently from other sites) we move 
min{X(x), K — X(x + 1)} particles with the largest indices from the site x to the site x + 1. 

Lemma 2.1 The order function I is preserved under the action of the map T. 

Proof. Straightforward. | 

For a given configuration X, associated to the collection of particles 5, for each particle 
£ G H we introduce the notion of velocity which is equal to 1 if the particle moves after 
the application of the map T or otherwise. Accordingly we shall say that the particle £ is 
free if v(£) > and jammed otherwise. Summing up velocities of individual particles we obtain 
moments (the total velocity of particles at a given site) of lattice sites in the configuration X: 

v(X,x) := "(£)■ 
\t\=x 

An immediate calculation shows that 



Lemma 2.2 v(X,x) = mm{X(x),X*(x + 1)}. 



From the point of view of the description in terms of individual particles we introduce the 
notion of the dual configuration X*(x) := K — X(x) for any x G 2Z 1 , which describes empty 
places in the original configurations of particles X. Therefore to describe the dynamics of empty 
places we consider the dual map T* : X^f — > Xq^ = (Xq")* whose action is defined by the relation 
T*X = (TX*)* and can be written as follows. 

Lemma 2.3 For X G Xq" we have 

T*X(x) = X(x) - min{X(x), X*(x - 1)} + mm{X(x + 1), X*(x)}. 

Observe that the dynamics of empty places is exactly the same as the dynamics of particles, 



but occurs in the opposite direction. Obviously both the above formula and the relation (2.1) 
describe the mass conservation rule: the number of particles at a given site in the new config- 
uration is equal to the number of particles at the same site in the original configuration minus 
the number of particles leaving it and plus the number of particles coming to this site. 

By a jammed cluster (of particles) we shall mean a locally maximal group of consecutive sites 
on the lattice containing at least one jammed particle at each site. Accordingly the jammed 
cluster in the dual configuration defines the cluster of free empty places in the original configura- 
tion. The locally maximal property means that any enlarging of the considered group breaks the 
definition, i.e. both immediate neighboring sites to the cluster do not contain jammed particles. 

Consider two subspaces of the space of configurations Xq' which shall play an important role 
in our analysis. The first of them is the space of configurations of free particles: 

Free(K) := {X G Xf : v(£) = 1 V£ G X}, 

and the second one is the space of (space) n-periodic configurations: 

Per(n, K) := {leX?: X(x) = X(x + n) Vx G 2Z 1 .} 

A trivial calculation show that both of these spaces are invariant with respect to the dynamics. 

Lemma 2.4 T : Fiee{K) -> Free(iC), T : Per(n,K) -> Per(n,K) for any n,K G 7L\. 



Proof. Observe that the restriction of the map T to the space of configurations of free 
particles is equivalent to the shift operator in this space, from where the first statement follows 
immediately. To prove the second statement notice that according to the formula (dynamics) 

TX(x) = X(x) + mm{X(x - 1), K - X(x)} - mm{X(x),K - X(x + 1)} 

= X(x + n)+mm{X(x-l + n),K - X(x + n)}-mm{X(x + n),K -X(x + l+n)} = TX(x + n) 
due to the n-periodicity of the configuration X G Per (n, K). | 



Denote by X = (a) the n-periodic configuration X G Per(n, K) consisting of the periodically 
repeating word a = a±a2 ■ ■ ■ a n with G {0, 1, . . . , K} and such that X{\) = a%, for example 
X = (1234) = . . . 123412341234 . . .. 

It is worth notice that despite the statement of the previous Lemma the minimal period 
of the configuration may not be preserved under the dynamics. Indeed, consider a 4-periodic 
configuration (1100) and observe that for K = 1 we have T(1100} = (1010) G Per(2, 1). 

To deal with more general and still statistically homogeneous configurations in the next 
section we introduce a more interesting subset of configurations - regular configurations for 
which as we shall show the statistical description makes sense. 



3 Space of regular configurations 



For a configuration X £ Xq- we define the notion of a subconfiguration X% := {X(k),X(k + 
1), . . . ,X(n)}, i.e. a collection of entries of X between the pair of given indices k < n, and 
introduce the corresponding density and the average velocity: 



where m(X^) := J2x=k-^( x ) stays for the number of particles in the subconfiguration X%. [[ 

By the density and the average velocity (of particles) of a entire configuration X £ we 
mean the following limits (if they are well defined) : 

p(X) := Jim p{X n _ n ), V(X) := lim\(X\), 

otherwise one can consider the corresponding upper and lower limits, which we denote by p±{X) 
and V±(X). 

Notice that both these quantities are well defined in the case of space periodic configurations 
(belonging to Per(i^)) in distinction even to the simplest case when a configuration X consists of 
free particles (i.e. belongs to Fiee(K)). Thus in the general case important statistical quantities 
p(X), V(X) may be not well defined. To be able to deal with the space of configurations satisfying 
a reasonable statistical description we introduce the following space of configurations. 

We shall say that a configuration X satisfies the regularity assumption (or simply regular) if 
there exists a number p £ [0, K] and a strictly monotone function ^(n) — » as n — > oo (which 
we call rate function), such that for any n G 5Z 1 , N £ 7L\ and any subconfiguration X^f-^ of 
length N the number of particles in this subconfiguration m(X^ti ) satisfies the inequality 



N 



P 



< <p(N). (3.1) 



It is clear that at for a configuration X satisfying the regularity assumption the density of 
particles p(X) is well defined and is equal to the value p in the formulation of the assumption. 
The space of configurations from Xq satisfying the regular assumption with the density p and 
the rate function ip we shall denote by Reg(p, cp, K). 

The main result of the paper formulated below describes the restriction of the dynamics to 
the space of regular configurations and will be proven in the rest of this section and the next 
one. 

Theorem 3.1 Let the initial configuration X £ Reg(p,(p) with p ^ K/2. Then after a finite 
number of iterations t < t c = t c (p,Lp) := ^(f' 1 ^ — P) + I) 2 f or the configuration T l X the 
average velocity of particles becomes well defined and is equal to min{l, ^ — 1}. Moreover for 
any t>t c we have T l X G Free(K) if p < K/2 and (T'X)* G Free(K) if p > K/2. 

To analyze properties of regular configurations we introduce the binary relation 'domination', 
which we denote by h, on the set of configurations Xq^ as follows: X h Y if and only if for any 
n £ 2Zi, N G 7L + there exists a pair n_, n+ G 7L such that 



1 Observe that in the definition of the average velocity we consider only particles from sites till n — 1. This is 
related to the fact that velocities of particles in the site n is not defined by the subconfiguration X^. 



Lemma 3.1 The relation h is an order relation, i.e. it is reflexive and transitive, but this 
relation is not symmetric. 



Proof. The proof of the reflexivity, i.e. that X \- X for any X G Xq" is straightforward. To 
prove the second statement consider a pair of configurations X h Y h Z. By definition for any 
n,k G 7L, N G 7L + we have 

Thus for any k, N there exists n_,n + ,n'_,n' + G TL such that 

{KlXr) < < m{Ztt») < m(Y<^ N ) < m(X^). 



m 



Therefore X h Z. It remains to prove the absence of symmetry, i.e. that there exists a pair 
of configurations X h Y such that the relation Y \- X does not hold. Let X(l) = 1, while 
= for all i / 1, and let Y(x) = for all x G 2Z. Then clearly Ih7 but the opposite 
relation does not hold. ■ 



Lemma 3.2 Let X \- Y and let X G Reg(/?, </?, if). Then Y G Reg(p, </?, K) as well. 

Proof. According to the definition for any n, N there exists a pair n-, n+ such that 

m(X: : +f ) < m(y£f) < ). 

Thus 

(m < mix?-:?) m {Y^ N ) < 
-V(N) < — p< j± p< ^ p < <p(N), 

which yields the desired statement. | 

Now we are ready to prove that the set of regular configurations is invariant under the 
dynamics. 

Lemma 3.3 T : Reg(p, tp, K) — > Reg(p, <£, if) /or any triple (p, ip, K). 

Proof. Let a configuration X G Reg(p, ip, K). We need to show that the configuration TX also 



satisfies the same assumption. To do it we shall prove that X h TX, from where by Lemma ^2 
we shall get the desired statement. Fix arbitrary integers neS and iV G S + and consider the 
sub configuration (TX)^t^ . The number of particles in this subconfiguration differs from the 
number of particles in the subconfiguration X™^± by the number of particles P coming from 
the site n to the site n + 1 and the number of particles P+ coming from the site n + N to the 
site n + N + 1, i.e. 

m((TX)£+f ) = m(X#f) + P_ - P+. 

There might be four possible situations: 

(a) X(n) + X(n + 1) < K and X(n + N) + X(n + JV + 1) < if. Then P_ = X(n), 
P + = X(n + N), and thus 

m((TX)^f ) = miXZl?) + X(n) - X(n + N) = m{X^ N - 1 ). 



(b) X{n) + X(n + 1) < K and X(n + N) + X(n + N + 1) > if. Then P_ = X(n), 
P+=K - X(n + N+1), and 

m(TOf) = m (Kti) + + * " X(n + N + 1) 
= m{X^ N ~ l ) + X(n + TV) + X(n + N + 1) - K > m{X^ +N ~ l ). 
On the other hand, 

m{{TXy$») = m(X^) + X(n) + K-X(n + N + 1) 

= m(^+f + 1 ) + X(n) +X(n + 1)-K< m(X£? +1 ). 

(c) X(n) + X(n + 1) > K and X(n + N) + X(n + iV + 1) < K. Then P_ = K — X(n + 1), 
P+ =X(n + N), and 

«n(Wf ) = m(X^) + K - X(n + 1) — X(n + N) 

= m{Xl +N ^) - X(n) +X(n + N)+K - X(n + 1) - X[n + N) > m{Xl +N - 1 ). 
On the other hand, 

m{{TX)%») = m(X^) + K - X(n + 1) — X(n + N) 

= m{X^ +1 ) -X{n + N+l) + K-X{n + N)< m(X^ +1 ). 

(d) X(n) + X(n + 1) > K and X(n + N) + X(n + TV + 1) > K. Then P_ = K — X(n + 1), 
P + = K - X(n + N + 1) and 

= ™PC+f) + K - X(n + 1) - K + X(n + N + 1) = m(^ +1 ). 

Therefore in all four possible cases we have found subconfigurations in X approximating (by 
the number of particles) those in TX from both hands, which yields the statement of Lemma. | 



Lemma 3.4 (Reg(p, tp,K))* = Reg(K — p, ip, K) . 

Proof. Let X G Reg(p, tp, K). Then for any n G 2Z 1 , N G ffi^ we have 

m((X*) Tt +f ) = m(«if) - X)-+f ) = K ■ N — m(X^). 



Therefore 



P ((x*r n x?)-(K-P) 



m((XT n ^ 
N 



(K-p) 



N 



P) 



< <P(N). 



Consider now the connection between spaces of periodic configurations and regular ones. 
Clearly, for any configuration X € Pei(n,K) the notion of density p(X) is well defined and 
p(X) = m(Xi )/n. To specify the density we denote by Per p (n, K) the set of configurations 
from Per(n, K) having the same density p. 



Lemma 3.5 For any p, n, K we have T : ~Pei p (n, K) — > Per p (n, K) and Per p (n, K) C Reg(p, p(l- 

JL\IL 

K> TV' 



Proof. For a given configuration X € Per(n, K) denote p := p(X) = m(Xf)/n. The first 
statement immediately follows from Lemma |2.4| and the fact that the number of particles on the 
period of the configuration cannot change under dynamics. Now each positive integer N can be 
represented as N = kn + I, where k € {0, 1, . . .}, I € 0, 1, . . . , n — 1. For any I < n — -j^n the 
number of particles in the subconfiguration X^f\ can be estimated from below as 

m{X%» ) > pkn. 

Therefore 

Otherwise, if Z > n — -jj^n we have 

m(Xl~\.\) > pkn + K(l -n + j^n) = pkn + Kl - Kn + pre. 

Thus 

m{X%Xx) pkn + Kl - Kn + pn 

9 N ~ 9 kn + l 



l(n-l)(K-p)<^(K-p) = ( p(N). 



Now we shall use estimates for the number of particles in the subconfiguration X^f\ from 
above. If / < -&n then 

m{X x x Xi) < pkn + Kl 

and 

m(Xg?) pkn + Kl 1 sP ( A n 



Otherwise, if I > fen then 



N r ~ kn + l r N x ' ' ~ N K r 'K 

K 



m{XlXi) < pkn + pn 
and 

HKXi ) / P^ + pn 1 , 

N 9 ~Tn~+T ~ 9 = N 9 ^ ~ ^ 

< Jl 9{n ~ K U) = ^ Pn(1 ~K ) = ^ N) - 



One can easily check that for any (space) periodic configuration the notion of the average 
velocity is well defined. It is of interest that for more general class of regular configurations this 
is not the case even for K = 1. Denote a = 1100 and b = 1010 and consider the configuration 
X constructed as follows: 

. . . bbbbaaaa bbaa ba ab aabb aaaabbbb . . . , (3.2) 

i.e. Xf = ab, X°_ 7 = ba, X| 4 = aabb, XZ23 = bbaa, etc. Notice that in each subsequent series 
the number of consequent elements aa . . .a and bb ... b doubles. 



Lemma 3.6 The configuration X defined as (3.&) is regular (X £ Reg(l/2, 1/N, K)), but the 
average velocity is not well defined. 



Proof. Observe that m(XHf k ) = 2k for any i,k, while 2k < m(xi+* k+j ) < 2k + 2 for any 
j G {1,2,3}. Therefore the configuration ( |3.2| ) is regular with the density 1/2 and the rate 
function tp(N) = 1/N. 

Let us calculate now the average velocity on various subconfigurations. First consider a sub- 
configuration starting from the 1st element and containing the full series aa . . . bb, i.e. X 1 
This sub configuration for any k contains the same number of elements a and b and hence 

v(x? 2k+1 -V) = 1.1 + 1.1 = -. 

y 1 ' 2 2 2 4 

Similarly due to the symmetry of the configuration X we have V \X°_ 2( ^ k+1 = 3/4 and 

therefore |V(X2" 2 t 2 +2 ) - 3/4| < 5 • 2~( fc + 3 ). Thus 

V(xf 2 t+ 2 ) -> 3/4 as k -> oo. 

Another type of subconfigurations that we consider differs from the previous one by the fact 
that it contains an additional (full) series of elements a ... a in the end, i.e. X 1 = 

^3-2 fc+1 -2 



V(Af 



1 . 2 fc+i + 3 . 2 ( 2 fc+i _ -q 2 . 2 fc+i _ 3 2 



y 2 fc + x + 2(2 k + 1 - 1) 3 • 2 fc + x - 2 3 

as k — > oo. Therefore using again the symmetry of the configuration X we have 

VprifJJii) 2/3 as fc^oo, 
and thus different subsequences of /c lead to different average velocities. 



4 Traffic jams and simple properties of the dynamics 

Recall that in a configuration X sites between x' and x" belong to the jammed cluster if for 
any integer x £ [x',x"] the inequality X{x) + X(x + 1) > K holds true. Similarly consecutive 
sites for which this inequality does not hold belong to a free cluster. The site x is called free if 
X{x)+X{x + l) < K. 

Lemma 4.1 For each configuration X E and for each site x € 2Z 1 we have 

min{X(x - l),X(x),X(x + 1)} < (TX)(x) < m&x{X(x - l),X(x),X(x + 1)} 
and thus for any t £ 2Z^ 

mm{X(x)} < mm{(T t X)(x)} < max{(T*X)(2;)} < max{X(x)}. 

Moreover, the upper and lower limits may be not preserved under dynamics: 

- 3X € Xq- such that max x {(TX)(x)} < m&x x {X(x)}; 

- 3X € Xq^ such that mm x {(TX)(x)} > m.in x {X(x)}; 

Proof. Fix a configuration X S Xj^ and a site x € 5Z 1 . Then denoting by P(x) := min{X(x — 
1),K — X(x)} the number of particles moving to the site x (from the site x — 1)) under the 
action of the map T we get 

(TX)(x) = X(x) + P(x) - P{x + 1). 



Consider all 4 possibilities: 

(a) X(x - 1) < K - X(x) and X(x) < K - X{x + 1). Thus 

(TX)(x) = X(x) + X(x - 1) - X(x) = X(x - 1). 

(b) X(x - 1) > K - X{x) and X(x) <K- X(x + 1). Thus 

{TX)(x) = X(x) + (K - X{x)) - X(x) = K — X(x) < X(x - 1). 
On the other hand, in this case 

(TX)(x) = K — X{x) > X(x + 1). 

(c) X(x - 1) < K - X(x) and X(x) > K - X(x + 1). Thus 

(TX)(x) = X(x) + X(x - 1) - {K - X(x + 1)) = X(x-l)-{K-X(x)) + X(x + l) <X(x+l). 
On the other hand, in this case 

{TX)(x) = X(x - 1) + X(x) — (K — X(x + 1)) > X(x - 1) + X(x) - X(x) = X(x - 1). 

(d) X(x - 1) > K - X(x) and X(x) > K - X(x + 1). Thus 

(TX)(x) = X(x) + (K - X(x)) — (K — X(x + 1)) = X(x + 1). 

Thus the first statement of Lemma holds true in all situations. 

It remains to construct examples of configurations satisfying the last two statements of 
Lemma. Let K = 2. Then T : (221022) -► (211122) and T : (002100) -► (001110). In the 
first example the minimal value becomes 1, while in the second example the maximal value 2 
becomes 1 under the action of the dynamics. | 

Introduce a map marking global maxima of a configuration M : ~Kq — > Xq as follows: 
MX(x) := 1 if X{x) = m.aXy{X(y)} and MX(x) := otherwise. We define also arithmetic 
operations with configurations X, Y G X^: 

(X + Y){x) := mm{X(x) + Y(x), K}, (X - Y)(x) := max{X(x) - Y(x), 0}. 

Using this notation we can formulate the following decomposition result. 

Lemma 4.2 For a given X £ Xf ifVx £ 7l\ such that MX (x) = 1 holds X(x-1)+X(x) < K 
and X(x) + X(x + 1) < K then 

TX = T(X - MX) + T(MX), 

otherwise 

TX = T ]k _ 1 (X - MX) +T h MX, 

where means the restriction of the map T to X^" -1 . On the other hand, there exists a 

configuration X G X^" such that 

{T*X)(x) + (Tf K1 X)(x). 

even if m&x x {X(x)} < K. 



Proof. The statement about the decomposition follows immediately from the definition of 
the dynamics while the following example demonstrate the second statement: XL : (1221) — ► 
(1212), Tj 2 : (1221) -» (2211). ' | 

These results demonstrate rather counter intuitive properties of the considered model of 
traffic flows. For example, from Lemma |4.1| it follows that if for a given initial configuration one 
traffic lane is not occupied (along the entire lattice), then these property holds for any moment 
of time. So it looks that the dynamics will not change if the road will be made narrower by one 
lane. However this is completely wrong, which was demontstrated in the second statement of 



Lemma 4.2 



Another example gives the following seemingly evident (but wrong) decomposition, which 



one would expect instead of the more complex decomposition described in Lemma 4.2. Assume 
that for a configuration X we have X(x) > 1 for all x £ 7L . Then it looks reasonable that the 
dynamics of the configuration, restricted to the lanes 2, 3, . . . , K should be the same as in the 
original one, i.e. 

T lK X = T lK _ 1 (X -(!}) + T lK ((l))- 
The following example of a periodic configuration shows that this is not the case: 

Tj 2 (1221> = (2211), T] 1 (0110) +Tj 2 (llll) = (0101) + (1111) = (1212). 

Lemma 4.3 Let X^f™ be a jammed cluster of length n in the configuration X. Then 

(TX)(x) = X{x + l) Vx e {k + 2,...,k + n}, (4.1) 

(TX)(k + n+l) = K-X(k + n + l), (4.2) 
(TX)(k + 1) = X(k) + X(k + l)+X(k + 2)-K (4.3) 
and if the site k — 1 does not belong to another jammed cluster, then 

(TX)(k-l) = X(k-2), (TX)(k)=X(k-l), (4.4) 

otherwise 

(TX)(k-l)=X(k), (TX)(k) = K — X(k). (4.5) 

Proof. First let us show that TX(x) = X(x + 1) for all x E {k + 2, . . . , k + n}. Observe that 
by the definition of a jammed cluster we have 

X(x-1) + X(x)> K, X(x) +X(x + 1) > K. 

Thus 

X(x -1) > K - X(x), X{x) > K-X(x + 1). 

Therefore after the application of the map T exactly K — X(x) particles comes to the site x 
from the site x — 1, while K — X{x + 1) particles leaves it. Therefore 

TX(x) = K- X{x) + X{x) -{K- X(x + 1)) = X(x + 1), 

which proves the equality (^^). Observe that this equality makes sense only if n > 2. 

The site k + n + 1 is the first free site after the jammed cluster. Therefore all particles from 
it moves to the site k + n + 2 under the action of the map T, while exactly K — X{k + n + 1) 



particles moves to the site k + n + 1 from the last site of the considered jammed cluster. This 
gives the equality (|4.2| ). Notice that from this inequality we get that 

(TX)(k + n) + (TX)(k + n + l) = K, 

i.e. the site k + n does not belong to a jammed cluster in the configuration TX. 

Clearly the site k cannot belong to another jammed cluster, otherwise the site k + 1 would 
not be the first site of the considered jammed cluster. 

By definition under the action of the map T all particles from the site k moves to the site 
k + 1, from where exactly X(k + 2) — K particles moves to the site k + 2. Thus we get the 
equality (4.3). 

Consider now the case when the site k — 1 does not belong to another jammed cluster, i.e. 
X(k — 1) + X{k) < K. This immediately gives the formulae ( |4.4| ) for the number of particles in 
the sites k — 1 and k. 

If n = 1 then 

(TX)(k+l) + (TX)(k+2) = X(k)+X(k+l)+X(k+2)-K+K-X(k+2) = X(k)+X(k+1) < K. 
If n > 2 then 

(TX)(k + i) + (TX)(k + i + 1) = X(k + i - 1) + X(k + i) 
for all i G {1, . . . ,n - 2}. Thus 

(TX)(k + n-l) + {TX)(k + n) = X(k + n) + X(k + n + 1) > K, 

(TX)(k + n) + (TX)(k + n + 1) = X(k + n+l) + K - X(k + n + l) = K. 

Therefore in both cases the site which was the last one in the jammed cluster becomes, 
and if n > 2 the site k + n — 1 turns out to be the last site in the jammed cluster in TX. 

Summarizing, in the case when the site — 1 is free we get the following representation for 
TX in the neighborhood of the considered jammed cluster: 

k-1 k k+1 k+2 k+3 . . . k+n-1 k+n 

k+n+1 

X= . . . X(k-l) X(k) [X(k+D X(k+2) X(k+3) ... X(k+n-l) X(k+n)] 
X (k+n+1) 

TX= ... X(k-2) X(k-l) (TX)(k+l) X(k+3) X(k+4) ... X(k+n)] X(k+n+l) 
K-X (k+n+1) ... 

By square brackets we denote the boundaries of jammed clusters. Observe that the first site 
of the new cluster might be either k or k + 1. 

In the alternative case when the site k — 1 is the last site of the previous jammed cluster, 
the representation differs only at sites k — 1 and k: 

k-1 k k+1 k+2 k+3 . . . k+n-1 k+n 

k+n+1 

X= . . . X(k-l)] X(k) [X(k+D X(k+2) X(k+3) ... X(k+n-l) X(k+n)] 
X (k+n+1) 

TX= ... X(k) [K-X(k) (TX)(k+l) X(k+3) X(k+4) ... X(k+n)] X(k+n+l) 
K-X (k+n+1) ... 

Indeed, applying the first statement of Lemma (which has been already proven) to the 
previous cluster, we get that (TX)(k — 1) = X(k). On the other hand, the number of particles 
moving from the site k — 1 to the site k is equal to K — X(k), while all the particles that were 
at site k move to the site k + 1 (since the site k does not belong to a jammed cluster). 



(TX)(k) + {TX){k + 1) = K- X(k) + X(k) + X(k + 1)+ X{k + 2) - K 



= X(k + 1) + X(k + 2) > K, 

therefore the new jammed cluster has the same length and is located at the sites from k to 
k + n — 1. ■ 

Corollary 4.4 For each jammed cluster (TX)^™ +1 of length n > 1 we have (TX)(k + n + 2) + 
(TX){k + n + 3) = K. 



Proof. Immediately follows from the equality (4.2) 



This implies that the distance between two consecitive jammed clusters is at least 2. 

Lemma 4.5 Let X^t^ be a jammed cluster of length n in the configuration X. Then neither 
its length, nor the number of particles in it cannot increase under dynamics. Moreover, if 
X(k — 1) + X(k) < K then number of particles in the jammed cluster decreases at least by 
K — (X(k — 1) + X{k)) > after the application of the map. 

Proof. Consider first the case when the site k — 1 does not belong to another jammed cluster, 
i.e. 

X(k - 1) + X(k) < K, X(k) + X(k + 1) < K. 

Clearly, in this case the site k—1 cannot be the first site of the jammed cluster in the configuration 
TX. Therefore in the worst case the jammed cluster is located at sites from k to k + n — 1, i.e. 



its length is at least not larger than of the considered one. Applying Lemma [L3| we can estimate 
from above the difference between the number of particles in the new cluster and the old one as 
follows: 

[(TX)(k) + (TX)(k + 1)} - {X(k + 1) + X(k + 2)] 

= X(k - 1) + X(k) + X(k + 1) + X(k + 2) - K - X(k + 1) - X(k + 2) 

= X{k - 1) + X(k) - K < 0. 

Hence the number of particles in this case cannot increase, and moreover this number decreases 
if X(k - l)+X(k) < K. 

It remains to consider the case when the considered jammed cluster is immediately preceded 
by another jammed cluster. Again by Lemma I 



(TX)(k - 1) + (TX)(k) = X(k) +K- X(k) = K. 
Hence the site k — 1 does not belong to the jammed cluster. On the other hand, 

{TX)(k) + (TX)(k + l) = K- X(k) + X(k) + X(k + 1) + X(k + 2) - K 

= X(k + 1) +X{k + 2) > K, 

since the site k + 1 belongs to the jammed cluster. Thus the site k + 1 is the first site of the 
jammed cluster in the configuration TX, which lies in sites from k to k + n, i.e. its length is 
exactly the same as of the old one. Applying the same trick as above to calculate the difference 
between the number of particles in the new cluster and the old one we get: 

{{TX)(k) + (TX)(k + 1)} - {X(k + 1) + X(k + 2)} 



= K — X(k) + X(k) + X(k + 1) + X(k + 2) - K - X(k + 1) - X(k + 2) = 0. 
Therefore even the number of particles in the jammed cluster is preserved in this case. | 
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Figure 1: Two examples of the dynamics of n-periodic configurations: (a) K = 2, n = 5, 
p = l = if/2, (b) K = 3, n = 5, p = 7/5 < if/2 = 3/2. 

Lemma 4.6 Let n € 2Z^ ; / o(X^ rt+1 ) < if/2, and let in the subconfiguration X^ n+1 there 
is at least one jammed site. Then there is an integer i £ {1, . . . , 2n — 1} such that X(k — i) + 
X(k + i + l) < K. 

Proof. Assume that this statement does not hold true. Then for any two consecutive sites x 
and x + 1 in this subconfiguration we have X(x) + X(x + 1) > if . On the other hand, for the 
jammed site y we get X(x) + X(x + 1) > if + 1. Thus 

In 

^X(k + x) > n + 1, 

x=l 

which contradicts to the fact that the density is less or equal to 1/2. | 

These results yield the following property: For any given subconfiguration the number of 
particles in any jammed cluster completely contained in this subconfiguration is a nonincreasing 
function of time and achieves its lowest possible level under dynamics. 

Lemma 4.7 Let X € Reg(p, ip, if) with the density p < if/2. Then after at most t c = t c (p, </?) = 
|(</? _1 (y — p) + I) 2 iterations all particles in T l X for t > t c will become free. 

Proof. According to the definition of regular configurations M(n) - the maximal number of 
particles in subconfigurations of length n of the configuration X for each n € 7L\ satisfies the 
inequality 

M(n) , s 

n 

Thus for any n > N c := ip~ l {^- — p) it follows that M{n) < n/(2K). By Lemma |4.6| in each 
subconfiguration of length n there is a pair of consequent sites whose total number of particles 
Q is strictly less than if. Consider the dynamics of this pair of sites. According to our previous 
results, while the site ahead of them is free these two sites will simply move one position forward. 



In the opposite case, when the next site is the first site of some jammed cluster by Lemma 
the number of particles in this cluster will decrease by if — Q > 0. On the other hand, free 
particles and jammed clusters move in opposite directions each with the velocity 1. Thus the 
maximum time between the consecutive meetings of a jammed cluster and a pair of consequent 
sites with the total number of particles less than if does not exceed n/2. Let n be the smallest 
integer larger or equal to N c . Since after each such meeting the number of particles in the 
corresponding jammed cluster decreases at least by 1 and since the number of particles in this 
cluster is less or equal to M(n) we get the following upper estimate of the transient period: 

M(n) ■ n n 2 1 , i . . , 9 

tc < V < -r- < -Atp-\K/2 -p) + if. 



Proof of Theorem 3.1. After the preparation made in Lemmas 3~2|-4.7 we are able to finish 
the proof of our main result. Indeed, in the case of a regular configuration X £ Reg(p, <p, K) 



with the density p < K/2 by Lemma 4.7 for any integer t > t c the configuration T l X consists 
of only free particles. In the opposite case, when p > K/2, we consider the dual configuration 
X* £ Reg(-fT — p, ip, K) (by Lemma [O]) and the since the action of the dual map is equivalent 
to the main one but proceeds in the opposite direction we get that T l X* £ Free(-ff). 

It remains to prove the statement about the average velocity of the configuration T l X for 
each t > t c . Again we start from the case of low density p < K/2. Since the configuration T l X 
consists of free particles, velocity of each particle is equal to 1. Thus 

n(^ B J = l Vn£Z^, 

which both shows that the average velocity s well defined and that V(T*A) = 1. Now if p > K/2 
the dual configuration to the configuration Y := T l X again consists of free particles. Hence 
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as n — > oo since the density of the configuration Y* is equal to K 



Observe that in the proof of Theorem |3.1| we actually derived an estimate of the length of 
the transient period as t < t c = t c (p, (f) := i( ( / 9_1 (T — P) + 1) 2 which goes to infinity as p — > 1/2. 
This is the reason why Theorem [O] does not cover the boundary case p = 1/2, which we discuss 
below. 



Theorem 4.1 Let the initial configuration X £ Reg(-|-,'/?, K) and let x'{t) < x"{t) be positions 



of two fixed arbitrary particles at the arbitrary moment t. 
subconfiguration X^$ converges to 1 as t 



Then the average velocity of the 



oo. 



Proof. Denoting p := -|- and choosing a positive integer M we consider a configuration 
M X obtained from the configuration X £ Reg(p, 99, K) by the following operation: for each 
integer k we remove from the configuration X the closest from behind particle to the position 
kM. For a given positive integer M any integer N may be represented as N = kM + I with 
I £ {-M,-M + 1,...,M - 1,M} and k £ ZZ 1 . Then 



m( 



' M v n+kM+l\ 
A n+1 ) 



m{X, 



n+kM+l\ 
n+1 ) 



k 



and thus 



TO 



On the other hand, 
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N' 



Therefore X £ Reg(p — jj, f + 775 K) and according to Theorem 3.1 after a finite number 
of iterations t c the average velocity of the configuration T tc ( M X) becomes equal to 1 (since all 
the particles in this configurations are free). 
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Figure 2: Long transient periods (t c = n- 
(a) p = 2 = if/2, (b) p = 15/7 > K/2. 



• 1) of n-periodic configurations with K = 4 and n = 7: 



Making an opposite operation, namely inserting a particle to the configuration X to the 
closest from behind to kM empty position for each integer k, we obtain another regular configu- 
ration +M X € Reg(/9 + jj, (p + jj, K). Again by Theorem 3.1 after a finite number of iterations 
the average velocity of this configuration becomes equal to 



K 



1 = 1 + 



KM 



1 



as 



M 



OG. 



Thus both (arbitrary close as M — > oo) approximations X to the configuration X have 
after a finite number of iterations (depending on M) the average velocity deviating from 1 by 
0(1/M). It remains to show that the average velocity of a subconfiguration of the configuration 
X can be estimated from above and from below by those from above approximations. Let X 
and Y be two configurations such that X(x) < Y(x) for all x and let x'(t) < x"(t) be positions 
of two fixed particles in the configuration X at the arbitrary moment t. Denote by y'(t) < y"(t) 
positions of the same particles in the configuration Y. Then 

v(<l?)>v(<«) 



for any moment of time t. Indeed, additional particles in the configuration Y present only 
obstacles to the motion of other particles, thus making the average velocity slower (or at least 
not faster). | 



In the case of (space) n-periodic configurations numerical examples below demonstrate a 
much better estimate of the transient period: t c < n— 1, but it is rather unclear if it is possible 
to generalize this result for more general regular configurations. 

Moreover, it turns out that even the upper (attainable) estimate of the length of the transient 
period for an initial configuration from Per p (n, K) is not monotonous on the length of the period 
n and heavily depends on its parity. The above example demonstrate the estimate t c < n — 1 for 



odd values of n. Now following mainly ideas proposed in [15| we shall show that this estimate 
is rather different for even values of the period. 



Lemma 4.8 Let X G Per p (2n, K) for some n > 1. T/ien i/ie length of the transient period 
t c < n. 



Proof. We introduce an operator G mapping the space of configurations Xq into the space 
of two-side sequences of real numbers defined as follows: 



x-l 



GX(x) := X(i) -(x- l)K/2, 

where we set J2ido = Y^!i=-j f° r an y positive integer j. 
One can easily show that for any x £ 7L l we have 

X(x) = GX(x + 1) - GX{x) + K/2, 

and if additionally X G Per p (2n, K), then 

x+2n-l 

GX(x + 2n) - GI(i) = ^ (X(t) - if/2) = 2nif(p - 1/2) (4.6) 

i=x 

and for any iEK 1 

G(rX)(i) = max{GI(x - 1), GJ(x) - K/2, GX(x + 1)}. (4.7) 

This yields that 

G(T t X)(x) = max{max{GX(x - t), GX{x -t + 2),..., GX(x + t)}, 

m&x{GX(x - t + 1), GX(x - t + 3), . . . , GX(x + f - 1)} - K/2}. 

As usual we consider three possibilities. 

(a) Assume first that p < K/2. Then from the equality (4.6) we get that GX(x + 2n) < 
GX(x) for each x. Then for t > n we obtain 

G(T t X)(x) = max{max{GX(x - t),GX(x -t + 2),.. .}, 

max{GX(x - 1 + 1), GX(x - t + 3), . . .} - K/2}. 

Thus G{T t+1 X)(x) = G{T t X)(x - 1) and T <+1 A"(x) = T\x - 1). Substituting these equalities 
into (14.71) we obtain 



= max{0, G(T t X)(x) - G(T t X)(x - 1) - K/2, G(T t X)(x + 1) - G(T*X)(x - 1)} 

= max{0, T l X{x - 1) - K, T l X(x - 1) + T*X(x) - K}, 

from where T t X(x) < K — T l X{x + 1) for any t,x. Therefore for t > n all particles in the 
configuration T*X are free. 

(b) p = K/2. Since in this case GX(x + 2n) = GX(x) for all x we get that for t > n 

G(T t X)(x) = max{max{GX(2), GX(4), GX(2n)}, 

max{GX(l), GX(3), GX(2n - 1)} - K/2} 

if x — t is even and 

G(T t X)(x) = max{max{GX(l), GX(3), GX(2n - 1)}, 



max{GX(2), GX(4), GX(2n)} - K/2} 

otherwise. Therefore G(T t+1 X)(x) = G(T t X)(x ± 1) and thus T t+1 X(x) = T l X{x ± 1) for all 
x G 2Z 1 , which yields that T l X G Free(K). 

(c) p > K/2. This case follows from the argument applied in the case (a), since we can 
consider the dual configuration X*, for which the density of particles is less than K/2. | 



Lemma 4.9 Let X G Per p (2n + 1, K) for some n > 1. Then the length of the transient period 
t c <2n + 1. 

Proof. Any configuration X G Per p (2n + 1, K) belongs also to Per p (4n + 2, K). On the other 



hand the number 4n + 2 is even and thus by Lemma |4.q we get the desired estimate of the 
transient period. | 

It is of interest that in the space periodic case we can give a more detailed information about 
the dynamics in time. 

Proposition 4.10 For each configuration X G Per(n, K) and any integer t > t c = n the 
sequence {T t X(x)}t is n-periodic on t for each x G ZZ . 

Proof. This result follows from the fact that for each t the configuration T l X is n-periodic 



in space and by Lemmas |L8| and |L9| the length of the transient period t c < n. Thus for t >t c 
the configuration T l X consists either of free particles (if p(X) < K/2) or its dual satisfies this 
property. Therefore T t+n X{x) = T t X(x) for any x G 7L . Notice that this period in time might 
be not minimal (which can be as small as 2). | 

Observe that this construction heavily depends on the periodic space structure of the con- 
figurations, which rules out the generalization for a more general situation. 



5 On the chaoticity of the dynamics 

In the previous sections it was shown that for sufficiently large time the dynamics occurs either 
in Free(K) or in (Free(-fiT))*, i.e. the corresponding dual configurations belong to this space. 
Therefore to study asymptotic (as time goes to oo) properties of the dynamics we consider 
its restriction to the space of configurations of free particles (which contains the union of all 
attractors of the map T restricted to the set of regular configurations with the density less than 
K/2). The following result shows that this map is chaotic in the sense that its topological 
entropy (see definitions in ]i~2[| ) is positive. 

Theorem 5.1 h top (T, Pree(iT)) = In (^±^ + \ + > f or an V K G The re- 

mainder term R(K) above satisfies the inequality \R(K)\ < 2. 

Proof. All the particles in configurations on the largest (i.e. containing all others) attractor 
are free, i.e. X(x) + X(x + 1) < K. Thus the action of the map is equivalent to the right 
shift map with the upper triangular transition matrix (i.e. all elements in the first line are 1, 
all but the last are l's in the 2nd line, etc.). It is well known (see, for example, |l~2|]) that the 
logarithm of the largest eigenvalue of this matrix gives the topological entropy of the right shift 
map and thus the topological entropy ht op (T) of the traffic flow. Therefore the representation 



of the largest eigenvalue of the transition matrix which we shall give below finishes the proof. | 



To simplify the notation we denote N := K+l. Let A(N) = (a^) be the NxN left triangular 
matrix, i.e. dij = 1 for all i + j < N + 1 and a%j = otherwise. This is a nonnegative symmetric 
matrix, therefore its spectrum belongs to the real line and its largest eigenvalue A max (^4(iV)) is 
positive. 

Theorem 5.2 \ max (A(N)) = f N + ~ + where the remainder term \R\ = \R(N)\ < 2. 

Proof. For an integrable function / € L 2 consider the operator 

Lf(x) := f 1 X f(s) ds. 
Jo 

According to || eigenvalues of the operator L ordered by their moduli are equal to 

(_l)fe+i 



A& :- 



(k- 1/2)tt' 



while e(x) := cos(^x) is the eigenfunction corresponding to the leading eigenvalue Ai. For 
simplicity we shall use the notation A := A(N), A := Ai = ^, := e(k/N), and e = i. Now 
since the function e(x) is analytical, decreases monotonically, and its second derivative satisfies 
the inequality | d ^ - 1 < it follows that for each &; = l,2,...,iVwe have 

(Ae) k = I e(s) ds + —J2 (e(-) - e (_ )) + j±, (5.1) 

i=0 







where the remainder term \R±\ = \R\(N)\ < ^ < 1. Thus, introducing the operator Gf(x) : = 
/(0) — /(l — x), we rewrite the last equality as 

{Ae) k = Le(k/N) + ^-Ge(k/N) + i? l£ 2 . 

On this step our aim is to show that there exists a function g G C^O, 1] orthogonal to e (which 
means that J g ■ e = 0) such that the following relation holds true: 

(L + £ -G)(e + = (- + s-)(e + £ -g) + R 3 e 2 , (5.2) 

Z Z 7T TT Z 

where again \R$\ < 1. 

Since the operator L is symmetric the orthogonal complement to the function e is invariant 
with respect to L. Thus for some constant a and a bounded function h € C 1 independent on e 
and orthogonal to e we have 

Ge = ae + /i. 

Let us calculate the constant a. Since h is orthogonal to e, then multiplying the both hands of 
the previous equality by e and integrating (observe that Jq 1 e • h = 0) we get 



On the other hand, 



e(x) ■ (e(0) — e(l — x)) dx = a e 2 (x) dx. 

Jo 

f 1 , , 1 Z" 1 ,7T , , 2 

/ cos — x) ax = -, / cos —x ax = — , 
/n V 2 ' 2' ./n V 2 ; vr' 



Jo cos ^2 X ' COS ^ 2 ' ~ X ^ 2 Jo 



■ ^ 1 

sin - = — . 

X 7T 



Thus 

2 1 1 

(---) = n«- 

1 1 2 

from where a = -. 

7T 

Therefore for any function g E C° we have 

e e^ _ e, e„ e 2 „ /2 e2\ e„ _ s e 2 



(L + ~G)(e + - 5 ) = £e + -L 5 + -Ge + -Gg = f - + -- 1 • e + -(h + Lg) + 



Comparing this relation with the equality (|5.2| ) we come to the conclusion that h + Lg = —g 
or g := (L — ~)h. Notice that the right hand side of the last expression makes sense since h is 
orthogonal to e. Thus 



g(x) = {L- -y 1 fe(0) - e(l - x) - -e(x) 

7T \ 7T 



From the first two leading eigenvalues, we deduce that the norm of the operator L — — restricted 
to the orthogonal complement to the leading eigenfunction can be estimated from above by 



w- = o— . Therefore 

07T 07T 



Combining above estimates we deduce that there exist two vectors v , £ £ 1R such that 

Av = (-N + -)v + £, |£|<t«M, 

which yields the statement of Theorem by the following a'posteriori matrix perturbation argu- 
ment | jT^| . Let the equality Av = p,v + £ be satisfied for a symmetric matrix A, two vectors 

6 H n with ||£|| < e\\v\\ and a scalar p. Then the closest to p eigenvalue A of the matrix A 
satisfies the inequality |A — p\ < e. 

Indeed, for p = A the inequality becomes trivial, while otherwise 

\\v\\ < \\(a - /.i)- 1 !! . ip - »i)v\\ = ■ m\. 

\p — A| 

Thus \p - A| < |${ < e. | 

As we already mentioned this result corresponds to the steady states of our model for the 
case of initial regular configurations with 'low' traffic. Notice that in the opposite case of 'high' 
traffic p > K/2 each jammed configuration is in one-to-one correspondence with its dual one, 



for which p < K/2. Therefore the statement similar to Theorem 5T holds in this case as well, 
moreover 

h top (T,Fiee(K) U (Free(JQ)*) = h top (T, Free (IT)). 

Observe that for any positive integer K the topological entropy for the map T is strictly positive, 
which yields the chaoticity of the map. 



6 Statistics of typical configurations 



In this section we shall derive statistical information about typical configurations of particles. 
For configurations X € Free(-ftT) let us denote by S(n, K) the total number of different subcon- 
figurations of length n € TL 1 ^. 

Lemma 6.1 S{n,K) = X^ ax (A(K + 1)) + o(X^ x (A(K + 1))), where A(N) is N x N left 
triangular matrix. 

Proof. Denote by Si(n, K) the number of subconfigurations of length n consisting of only free 
particles and starting with the symbol i £ {0, 1, . . . , K}. Then we have the following recurrence 
relation: 

K-i 

Si(n + 1,K) = Sj{n,K). 

3=0 

The number of particles in each site of a configuration X may vary from to K, i.e. it may 
admit N := K + 1 different values, and the only additional relation that should be satisfied is 

X(x) + X(x + 1) < K Vx G {1,2, ... ,n}. 

Therefore these configurations are completely described by N x N left triangular transition 
matrix A = A(N). Thus we get S(n,K) = Yli=o Si(n, K), from where and from Theorem 5.2 
the statement of Lemma follows. ■ 

We shall say that a subconfiguration is blocking (non-blocking) if it contains (not contains) 
the symbol K. Then the number of non-blocking subconfigurations of length n is equal to 
S(n, K — 1). The fraction of blocking subconfigurations of length n is equal to 

S(n,N) - S(n,N-l) _ 1 _ S(n,N- 1) 
S(n,N) ~ S(n,N) 

Applying the asymptotic representation for the leading eigenvalue of the matrix ^4(iV) we get 
the following estimate: 

/ 2 1 1 \ " 

S(n,N)=(-N + - + o(-" 

\ 7T 7T 1\ 



Therefore 



S(n,N-l) _ / |(iV-l) + i+o(^) \ n _(,_]_, 
S(n,N) -{ fiV + i + o^) ) N +0[ N ) ) 



To derive a more deep information about the statistics and to be able to deal with periodic 
configurations one can use an approach based on Markov chain approximations. 

For for a given integer i £ {0, 1, . . . , K} and a configuration X £ Free(iT) denote by 7fj(X™) 
the fraction of sites x £ {1, 2, . . . , n} where X{x) = i, i.e. 

TT^Xf) := ~ #{x e {1, . . . , n} : X{x) = i}, 

while by 7r 4 - we denote the average of these fractions over all possible different subconfigurations 
of free particles of length n: 

_(n) _ SxgFree(X) 
' T; : Sill. I\) ' 



i 



Lemma 6.2 -k^ — > K ^}^ 1 -jrr^ as n — > oo. 



Proof. Continuing the same argument as in the proof of Lemma 6.1 we see that each config- 
uration can be considered as a realization of a Markov chain with N = K + 1 states numbered 
as 0,2,... , K and the following transition probabilities: 



Pij 



_ = {V{X(x + l)=j\X{x) = i} = ^X ri ]£j<K- 
1 otherwise. 



Clearly the iV-th power P of the transition matrix P = (pij) is strictly positive and thus the 
Markov chain is ergodic. Denote by m, i € {0, ...,K} its stationary probabilities, i.e. the 
probability to have X(x) = i. Then these quantities should satisfy the following system of 
equalities: 

K-i 

ETTj 



from where 

I + 1 

for all i = 0, 1, . . . , K. Solving the last system of difference equations we get 

K + l-i 

TTi = 7Tn, 

K + l u ' 

and eventually (since they sum up to 1) we come to 

_ K+l-i 2 
n ~ K + l K + 2 



To study statistics of n-periodic configurations of free particles one should take into account 
that there is an additional constraint: X(l)+X(n) < N. Denote by S(n, N) the total number of 
sub configurations of length n consisting of free particles and satisfying this constraint. Therefore 
the fraction of nonadmissible n-periodic configurations (which do not satisfy above constraint) 

K K K I K-i \ 

i=0 j=K-i i=0 \ j=0 / 

is asymptotically (for large K) equal to 

/ tt(x) [ vr(y) dy dx = \, 
Jo Jl-x o 

where tt{x) := 2(1 — x) corresponds to the limit (as K — > oo) distribution. 

Using these statistics one can easily obtain all correlation functions and to study large devi- 
ations. For example, we get that the fraction of blocking configurations in Free(if ) is equal to 
ttk = 2/{{K + 1){K + 2)), which in terms of traffic estimates how often the road with K lanes 
is completely blocked by moving cars. 
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